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<D 

-^-We summarize recent theoretical results for the signatures of strongly correlated ultra-cold fermions 
i n optical lattices. In particular, we focus on: collective mode calculations, where a sharp decrease 
^_dn collective mode frequency is predicted at the onset of the Mott metal-insulator transition; and 

(^correlation functions at finite temperature, where we employ a new exact technique that applies the 

^stochastic gauge technique with a Gaussian operator basis. 

T3 

qL Introduction 

o : 

Recent progress in obtaining bosonic and fermionic quantum gases has ush- 
^Hered in a new experimental paradigm in physics. An essential feature in these 
^^developments is that, although technically challenging, ultra-cold atom exper- 
^Jments have remarkably simple physical conditions. The significant issues are 
f"-~-that: 
■ 

<0» the underlying interactions are well understood, 

there are few parameters needed 
^» interactions can be tuned in many cases 

r^This will greatly help our understanding of many-body physics, since precise 
(theoretical models can be readily obtained at low density, without the compli- 
Ocations of a crystal lattice and reservoirs - as is so often the case at present. 



Thus, one can apply simple theoretical models with high accuracy, allow- 



o. . 

(jing innovative experimental tests of methods in quantum field theory (QFT). 
^^This could lead to novel fundamental experiments in physics, such as tests of 
: '"[massive particle quantum entanglement, as well as applications to metrology. 
r> \ An important element in recent investigations is the technique of coherent 
C^transformation of cold atoms to molecules, either using photo-association [1] 
or Feshbach [2] resonance methods. This allows a strong, tunable interaction 
between the atoms. Experiments have confirmed coherently coupled atom- 
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molecular models [3], with renormalization [4]. 

Bosonic molecular experiments include 133 Cs, 87 Rb, and 23 Na [5]. Fermi 
gases of 40 K and 6 Li atoms have proved to be ideal experimental systems ow- 
ing to their long molecular lifetimes [6]. These have resulted in cold molecule 
production [7], molecular Bose-Einstein condensation [8], and preliminary in- 
dications of fermion superfluid behaviour [9, 10]. 

The properties of ultracold Bose gases in optical lattices [11] are another 
'hot' topic, as a starting point for the exploration of strongly correlated many- 
body systems [12, 13]. Fermionic experiments in optical lattices are already 
underway [14], thus allowing the direct realization of the strongly coupled 
lattice of fermions known as the Hubbard model. In this paper we review 
some recently obtained results for the Hubbard model, in particular: 

• the use of collective modes as signatures for the metal-insulator phase- 
transition 

• a new exact Gaussian technique for calculating finite temperature correla- 
tion functions, without the Fermi sign problem of traditional Monte-Carlo 
methods. 



2 Hubbard model 

A one-dimensional, sinusoidal optical potential is easily created for trapped 
atoms. The results are different for fermions and for bosons, since because of 
the Pauli exclusion principle, at least two distinct spin eigenstates are needed 
to have on-site interactions with fermions. 
This leads to the (fermionic) Hubbard model, 



where Vj includes the external trap potential which we assume is harmonic, 
U gives rise to on-site interactions, and Uj describes tunneling between the 
sites. This model of interacting fermions is also used in high-T c superconduc- 
tivity, where it is more qualitative than quantitative in its applicability. The 
dimensionality and lattice structure of the model determine the tunneling ma- 
trix tij. In most of the paper, we will assume this is a one-dimensional lattice 
with = i<5j+ij, and Vj = mui^d 2 j 2 /2, where m is the fermion mass, luq the 
trap frequency, and d the well-separation. We also give correlation results for 
two-dimensional lattices in the last section. 




(1) 



j>=±i 
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Figure 1. Left graph: schematic illustration of effective single-particle band structure in a 
harmonically trapped fermionic Hubbard model. Right graph: phase-diagram vs the coupling (k ) 
and filling [y ). For details see Section 12.21 . 

2.1 Band structure 

In a single-band model as given above, one can expect a band insulator to 
form in the non-interacting limit, when all available sites have been occupied 
(i.e., two particles per site). However, this is complicated by the existence of 
the external potential. For a large, deep trap, the density increases towards 
the centre, so that the insulating region becomes localised, coexisting with a 
conducting region in the wings. 

When there are strong repulsive on-site interactions, the system develops 
a new band structure in which an insulating region forms at half- filling, i.e., 
with only one particle per site. This is shown in Figure (jT|). 

2.2 Local Density Approximation 

In the local density approximation (LDA), we assume that the trap is slowly 
varying compared to all correlation length scales. This allows us to use the 
piecewise exact solution for the uniform ID Hubbard model [15] to describe 
the strongly-correlated interacting ground state. 

It most useful to introduce an effective mass and frequency, corrected for 
the band-dispersion, as well as a dimensionless coupling constant and filling 
factor, where: 

• Effective mass, frequency: m* = h 2 / (2id 2 ), oj = Hujq (m/m*) 1 ^ 2 jt 

• Coupling constant k = U 2 / (8t 2 Nto) 

• Effective filling factor v = y2Nuj/n 



The effective filling factor can be understood as corresponding to a single par- 
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ticle per site (n = 1) at the trap centre for v = 1, although the interpretation 
is more complex in general, due to the effects of the trap potential. Note that 
v is defined as a global parameter for the entire trap, while the occupation per 
site (n) varies with radius. 

With these parameters, and using the LDA, with a fixed global chemical 
potential and harmonic trap potential, one obtains a phase-diagram [16, 17] 
as shown in the right-hand diagram in Fig (^Q). The labelled regions are as 
follows: 

A: a pure metallic phase. In this case there are no filled bands, and the quan- 
tum gas behaves analogously to a metallic conductor, with n < 2 everywhere. 

B: a single Mott insulator domain at the center, accompanied by two metallic 
wings. Here the strong particle-particle interactions prevent occupation num- 
bers larger than one, even though these are allowed by the Pauli principle. 

C: a metallic phase at the center surrounded with Mott insulator plateaus. 
In this larger particle number opens up a second band in the trap 

center, in which higher occupations with n > 1 are possible. 

D: a band insulator at the center with metallic wings outside. This is similar 
to phase B. Here, however, due to weaker inter-particle interactions there is 
no Mott phase, and the central band insulator with n = 2 is caused purely by 
the Pauli exclusion principle. 

E: a band insulator at the trap center surrounded by metallic regions, in turn 
surrounded by Mott insulators. In this case there is a large particle number 
and strong interactions, so both Mott insulator and band insulator regions 
can occur at different radii from the centre. 



3 Collective mode frequencies 
3.1 Luttinger approximation 

In order to describe the collective modes, we can make use of the Luttinger 
long- wavelength Hamiltonian [13], valid for small displacements. Because of 
the spin-density separation in the Hamiltonian, there are two types of waves 
present, described by the indices v = p, a. We assume here that the external 
potential only depends on the density, and does not transform one type of 
wave into the other, so that the inhomogeneous Hamiltonian is: 



Here the density and spin- wave velocities are u p , -Uo-respectively, and these 
have Luttinger exponents K p , K a . Using the LDA, we can therefore solve the 
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Figure 2. Breathing (ujb) and dipole (u>d) mode frequency vs filling factor v for zero, intermediate 

and strong couplings. 



resulting wave-equations for the collective mode frequency, with zero-current 
boundary conditions at the boundaries where the density vanishes [18]. 

The results are shown in Fig ©, in which we have solved for the collective 
fermionic modes in a trapped one-dimensional two-component Fermi gas, to- 
gether with an optical lattice [17]. These graphs show a sharp frequency dip as 
an unmistakable signature of Mott metal-insulator transition (MMIT) physics. 
This occurs just at the filling factor where the insulating region starts to form 
at the trap centre. 

The reason for this is simple to understand. Collective modes involve large 
density currents, which propagate through the lattice with a characteristic 
velocity u p . However, these velocities tend to zero in the neighborhood of an 
insulating region, and must be exactly zero inside the insulator. This effect 
slows down transport processes and therefore reduces all collective mode fre- 
quencies, which scale as velocity /length. The effect is strongest just when the 
insulating region starts to form, leading to a large fraction of the conduct- 
ing region having a low average velocity. As the insulating region grows, the 
part of the trap that is conducting is localized to the wings, so the frequency 
increases again due to the smaller characteristic lengths involved. 

However, there are some limitations to these methods. In particular, this is 
a linearized method, and hence is only valid for small displacements; it also 
only applies to zero temperature. 

Preliminary results in the non-interacting case using the Boltzmann equa- 
tions (circles in Fig (j2J) show that there are additional damping effects that 
are proportional to displacement. This leaves an unsolved problem, of how to 
deal with finite temperatures and large displacements in the interacting case. 
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4 Quantum simulation with Gaussian operators 

We now turn to the question of how to carry out many-body quantum field 
theory calculations without approximations. Quantum Monte-Carlo (QMC) 
techniques for fermionic problems are strongly limited by the Fermi sign prob- 
lem, which causes large sampling errors [19]. 

Here we give some results using a novel method constructed from a basis of 
Gaussian operators, which treats covariances (Green's functions) as phase- 
space variables. In principle, this method can simulate both fermions and 
bosons [20], including either thermal ensembles and dynamics. We will treat 
the case of finite-temperature fermionic ensembles for definiteness. 

4.1 Gaussian expansion 

Our approach is different from traditional_QMC methods. We expand the state 
density operator p in an operator basis A: 



where P( A ) is a probability distribution, sampled stochastically over the 
variables A which constitute a phase-space. Time (or inverse temperature) 
evolution is obtained by a mapping procedure which maps operator equations 
into stochastic equations: 



These equations are typically not unique, and have a 'stochastic gauge' 
degree of freedom, which must be chosen to make the distribution compact 
enough to eliminate boundary terms and minimize sampling errors [21]. 

4.2 Gaussian operators 

The mapping procedure depends on the operator basis chosen, and we shall 
choose a Gaussian operator basis, defined for both Fermi [22] and Bose [23] 
cases, as: 




(3) 



dp/dt = L[p] 



d A /dt = A + B C • 



(4) 



A±(n) = Q |I ± n| Tl : exp a 



(iTl-tlin]- 1 )^ 



(5) 



Here, annihilation and creation operators are included as a vector: a = 
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(Si, ...,(im) T ■ The upper signs apply for bosons, the lower signs for fermions. 

A typical application of these methods is to calculate thermal equilibrium 
correlation functions. Since the fermion number is usually unknown, the grand 
canonical distribution, p = exp(—(H — pN)r), is the most appropriate choice 
of ensemble. Here, p is the unnormalised density operator, r = l/ksT is the 
inverse temperature, and p is the chemical potential. The phase-space variable 
n is a complex matrix, which is best thought of as a (finite temperature) 

stochastic Green's function. Thus, /ajoj-y = (On^) p( T ), where {.) p indicates 

a phase-space average over the distribution P. 



5 Fermionic correlations 

Rewriting the canonical density operator as an equation for inverse tempera- 
ture evolution, one obtains: 



dp__ _1 



H-pN),p . (6) 



After applying the appropriate identities, one then obtains the following 
Stratonovich stochastic gauge equations for the Hubbard model, where O is a 
weight and n a is the stochastic Green's function for spin a = ±1: 

ar 

^ = I{(I-n CT )Ti 1 W + n CT Ti 2 )(I-n CT )} . (7) 

Here we have introduced a transition matrix T„ and real Gaussian noise 
terms defined by: 

• T-matrix: = % - 5 hJ {\U\{sn jtjj - a - n jtj , a + \) - p + ( j( s + 1 )/ 2 ^ r) } , 

• Noises: ((ff \t) ^V)) = 2\U\6(t - T')6 j;f 5 r>r , . 



5.1 Finite-temperature correlations 

The left graph of Fig @ shows the second-order correlation function 
for a simulation of the uniform one-dimensional Hubbard model, with both 
attractive and repulsive on-site interactions, compared to the analytic result 
at zero temperature [15]. 
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Figure 3. Hubbard models vs inverse temperature r. Left graph shows the second-order 
correlation function 92(0) = (rijj+njj—) / {njj^.) (fijj—) for a 100-site ID lattice at half-filling. 
\U\ = 2 and t = 1. The right graph shows a 16 X 16 2D lattice for chemical potentials /i = 2 (solid), 
fj, = 1 (dashed) and fi = (dot-dashed). Curves without crosses give the number of particles per 
site for fi = 1 (dashed) and fi = (dot-dashed). U = 4, t = 1, and 50 paths initially. 



The right graph shows the mean energies and particle numbers in a 2D 
lattice for a variety of chemical potentials, with no evidence for a Fermi sign 
problem that would cause enhanced sampling errors at low temperature. We 
note that this extension to higher dimensions is especially troublesome for 
traditional QMC techniques in the repulsive case. Our technique shows no ev- 
idence of increased sampling error at increased dimensionality, or with changes 
either to the filling factor or the sign of the coupling constant. 



6 Summary 

We have reviewed some recent theoretical developments for ultra-cold 
fermionic atoms in an optical lattice. 

For the case of the Mott metal-insulator (MMIT) transition in one dimen- 
sion, at zero temperature, we give a universal phase diagram based on the 
exact ID solutions, and valid in the LDA limit of large traps. A crucial exper- 
imental signature was predicted, namely a strong frequency dip at the filling 
factor corresponding to the onset of the MMIT at the trap centre. 

More generally, a novel exact technique for general dynamic and static 
Fermi calculations was presented. This can be used to calculate correlations 
of strongly-correlated fermions at any temperature - and for Hubbard models 
in any dimension. This is therefore a solution to the long-standing Fermi sign 
problem in the Hubbard model. 
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